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We study the electronic spin relaxation effect in the hole-doped monolayer and bilayer transition- 
metal dichalcogenides in the presence of the crystal defects. We consider realistic models of the 
lattice vacancy and actually estimate the spin relaxation rate using the multi-orbital tight-binding 
model. In the monolayer, the spin-relaxation time is found to be extremely long compared to the 
momentum relaxation time, and this is attributed to the fact that the spin hybridization in the band 
structure is suppressed by the mirror reflection symmetry. The bilayer TMD has a much shorter 
spin relaxation time in contrast, and this is attributed to stronger spin hybridization due to the 
absence of the mirror symmetry. 


PACS numbers: 72.25.Dc,73.63.Bd,85.35.Ds 

Monolayer transition-metal dichalcogenide(TMD) is a 
atomically thin two-dimensional semiconductor with a 
strong spin-orbit interaction^. In the electronic band 
structure, the spin and valley {K and K') degrees of free¬ 
dom are intercorrelated because of the broken inversion 
symmetry in the atomic configuration, implying various 
spin-dependent phenomena and potential applications to 
the spintronic devices^P— and electronicsiSrJ^. 

While the spin-orbit interaction is a key to control the 
electronic spin, it also causes the spin relaxation at the 
same time in presence of the impurity scattering. In 
the conventional semiconductors, the spin polarization 
of the conduction electron rapidly decays in various pro¬ 
cesses such as Dyakonov-Perel mechanism and Elliott- 
Yafet mechanismJi^— In the TMDs, on the other hand, 
the recent experiment showed that the carriers in the va¬ 
lence bands have a relatively long spin-relaxation time 
compared to conventional materials^!. The long spin life 
time comes from the peculiar band structure of TMD, 
where the spin-up and spin-down subbands split in oppo¬ 
site direction between K and K' valleys with a relatively 
large splitting width of the order of 100 meV [Fig. [1]. In 
the moderate hole concentration, therefore, the carriers 
at K and K' are fully polarized to opposite spin direc¬ 
tions, and then the inter-valley scattering is necessary for 
the spin relaxation. 

In this paper, we study the inter-valley spin relaxation 
effect caused by the non-magnetic short-range scatter- 
ers in hole-doped monolayer and bilayer TMDs. In the 
literature, the spin-relaxation in monolayer TMD was 
theoretically investigated for the conduction / valence 
bands and various different relaxation mechanisms 
The spin relaxation by the inter-valley scattering has also 
been discussed using effective impurity models^i^ while 
the spin-flip mechanism in actual crystal defect or impu¬ 
rity in the lattice structure has not been studied well. In 
what follows, we consider specific atomic defects^^ shown 
in Fig. [2] which are actually observed in real systems, 
and estimate the spin relaxation rate using the multi¬ 
orbital tight-binding model. There we show that the 
spin-polarized valence bands at K and K' are actually 


not pure spin states but includes small components of 
the opposite spin due to the spin-orbit interaction, and 
this enables the inter-valley scattering under the non¬ 
magnetic defects. In the monolayer TMD, such a spin 
hybridization is relatively weak because of the limitation 
by the mirror reflection symmetry ah, leading to a ex¬ 
tremely small spin relaxation rate (long spin relaxation 
time) compared to the momentum relaxation rate. In 
the bilayer TMD, on the other hand, the spins are more 
strongly hybridized due to the absence of the mirror re¬ 
flection symmetry, and as result, the spin relaxation rate 
is shown to be greater than in monolayer by the factor 
of 10^. 
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FIG. 1. Schematic description of the inter-valley spin scat¬ 
tering in monolayer TMD. 

In TMDs, the electric states around the Fermi energy 
consists of the d-orbitals {d^r'^-z^,dx'^-y 2 ,dxy,dxz,dyz) 
on transition-metal atoms (M = Mo, W) and the p- 
orbitals {px,Py,Pz) on chalcogen atoms (X = S, Se, 
Te). As the primitive unit cell contains a single 
transition-metal atom and two chalcogen atoms, we con¬ 
sider eleven atomic orbitals for each spin in a unit cell. 
The monolayer TMD has the mirror reflection sym¬ 
metry ah with respect to the layer plane, and there¬ 
fore the atomic orbitals can be classified into even or¬ 
bitals {d^z^-r^,dx^-y^,dxy,pt,Py ,Pl): ^-^d odd orbitals 
{dxz,dyz,Px ,Py ,pt), with respect to the eigenvalue in 
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FIG. 2. Atomic structure of (a) a single chalcogen-site va¬ 
cancy (V5c), and (b) multiple vacancies on a transition-metal 
sites and three chalcogen sites (Fmxs), in monolayer MX 2 . 


inverting the z-direction (out-of-plane direction). Here 
P/x is the superposition of the atomic orbitals 
on the vertically located pair of two chalcogen p-orbitals 
with and pj) for top and bottom. In the absence of the 
spin-orbit interaction, the even and odd orbitals indepen¬ 
dently form the even and odd energy bands, respectively. 
Fig. [3] shows the band structure of M 0 S 2 calculated by 
the first-principle method of quantum-ESSPRESSO^^. 



FIG. 3. Band structure of monolayer M 0 S 2 without spin- 
orbit interaction. Even and odd bands are indicated by green 
(dark-colored) and orange (light-colored) curves, respectively. 

To describe the motion of electrons, we derive the 
tight-binding model from the first principle band struc¬ 
ture by using WannierQO, which is a numerical package 
for minimized wannnier functions^. We first create the 
tight-binding Hamiltonian disregarding the spin-orbit in¬ 
teraction, and then introduce L ■ s spin-orbit term to the 
transition metal atoms. The Hamiltonian for the spin- 
independent part is written as 


cell i, respectively, and with s = ± represent the spin 
degree of freedom. The spin-orbit interaction for the d- 
orbitals of transition-metal atom is described in the basis 

(^32^— ^xy-! dxz-i ^yz^ a^^^" 
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Here the first three orbitals (daz^-r^, dxy) are 

even and the latter two {dxz, dyz) are odd under the 
mirror reflection ah. Since the spin flipping terms {ax and 
ay) always appear in the off diagonal block connecting 
even and odd orbitals, the spin orbit Hamiltonian can be 
written by 


Hso = X] <+ bfu°{s)b^, 

s,i 

+ -I-h.c , 


(3) 


where s represents the opposite spin to s, it® and u° rep¬ 
resent the diagonal blocks for even and odd states, re¬ 
spectively, and is the off-diagonal block. We neglect 
the spin-orbit interaction in the chalcogenide atom which 
is much smaller than that in the transition-metal atom. 

Before the detailed numerical calculation, we consider 
the mechanism of the spin-relaxation through the inter¬ 
valley scattering in the hole-doped TMD monolayer. It 
is similar to the conventional Elliot-Yafet process, while 
the even-odd classification due to the reflection symme¬ 
try imposes a limitation to possible scattering process. 
In TMD, the low-energy spectrum near the band gap 
is dominated by the even-orbital bands. The diagonal 
part of Hso is responsible for the band splitting between 
spin-up and spin-down branches due to icrz termsi. As 
schematically shown in Fig.[Tl the valence bands are spin- 
split in opposite directions between two valleys due to 
the time-reversal symmetry, and in the moderate hole- 
doped regime considered in the following, the Fermi en¬ 
ergy crosses only (iF, t) and (AT',4-) branches. 

The off-diagonal part of Hso hybridizes these low- 
energy even-state bands with the odd-state bands far 
from the Fermi energy with the opposite spin. The states 
near K and K' points are then expressed in the first-order 
perturbation as 

|A:,t) ~ I even, t) + ^|odd,;) 

\K',i) ~ I even, 4 .) -k ^|odd,t), (4) 


Ho = Y. + bf t°(*, j)b- (1) 

where a® and b® are the annihilation operators for the 
six even orbitals and the five odd orbitals in the unit 


where |even(odd), s) represents the eigenstate in the ab¬ 
sence of spin-orbit coupling, which is a direct product 
of the even (odd) orbital state and the pure spin state 
with s. Here A is the energy scale of the coupling ma¬ 
trix elements, and AE is the typical energy distance from 
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the Fermi energy to the odd state bands. We can show 
that A vanishes right at K and K' points^^, and linearly 
increases as the wave number is shifted from the valley 
center. Apart from the numerical factor (of the order of 
1 ), it is roughly written as 

A ~ \ka, (5) 

where A is the spin-orbit coupling constant in Eq. k 
is the relative wave-vector from K or K' (omitted in Eq. 
O), and a is the atomic scale constant. Fig. U shows 
the relative amplitude of spin-down component included 
in the up-spin state in M 0 S 2 monolayer (indicated by 
red dots), which is averaged over the Fermi surface and 
plotted against the corresponding hole density nh- We 
see the amplitude is roughly proportional to oc fc in 
consistent with Eq. ©• 


inter-valley spin relaxation rate (inverse spin relaxation 
time) is 

T,^27t\{K' imKt)fn,pS^ ( 8 ) 


where p is the density of states at the Fermi energy per 
valley and per area. 

On the other hand, the momentum scattering rate is 
dominated by the intra-valley scattering process without 
spin flip, and it is written as 

rp^27r|(iFt|Ff.|if t)|Vp52, (g) 


where {K f \Hi\K f) Va^jS. The ratio of two scat¬ 
tering rates then becomes 


Ts / A A' 2 A2 

-.^VKe) 


( 10 ) 
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FIG. 4. Relative amplitude of spin-down component in¬ 
cluded in the up-spin state in M 0 S 2 monolayer and bilayer, 
as function of the hole density Uh- 


We assume the scatterer Hi is non-magnetic and does 
not flip spins. Then the scattering matrix element be¬ 
tween |iF, t) and |A"',4-) is written as 

{K',i mK,t) = ^[(odd,t |A,|even,t ) 

-I-(even,! |i/i|odd,4, )]. (6) 

We see that Hi is required to break the mirror reflection 
symmetry to hybridize the odd and even orbital states. If 
we consider an atomic-scale scattering potential breaking 
mirror symmetry, we roughly have (even, s|FIi|odd, s) ^ 
Va^/S, where V is the energy scale of the potential and 
S is the total system area. As a result, the inter-valley 
matrix element can be estimated as 

{K' (7) 

When V is small and the scatterers are distributed 
sparsely in the system with a number density rii, the 


where we used A ~ Xkpa, and nt = 27rfc|./(27r)^ is the 
number of holes per unit area. The linear dependence 
of Fg/Fp on the hole-density Uh originates from the k- 
dependence of A. 

For M 0 S 2 , the strength of the spin-orbit interaction 
A is given by A « 0.073 eV— , and the odd-even energy 
distance AE is typically 2 eV. The lattice constant is 
a ~ 0.3 nm. At the moderate hole density of Uh = 10^^ 
cm“^, for example, the ratio Fg/Fp becomes 10“^. When 
assuming the disorder scattering rate Fp ~ O.OI eV, for 
example, we have the spin-relaxation time can be esti¬ 
mated to Tg = h/Ts 10 ns. 

The realistic short-ranged scatterers such as atomic de¬ 
fects or vacancies cannot be treated as a perturbational 
potential, and then we need to actually solve the scatter¬ 
ing problem to estimate the spin relaxation rate. We con¬ 
sider two specific mirror-symmetry-breaking configura¬ 
tions which are actually observed in the real materials^^: 
a single chalcogen-site vacancy (Vx), and a multiple va¬ 
cancy on a transition-metal sites and three chalcogen 
sites ( 141 X 3 ), which are shown in Fig. [5] (a) and (b), re¬ 
spectively. We consider a system with a single defect of 
each type in a tube geometry with a sufficiently large 
diameter, and calculate the transfer matrix T using the 
tight-binding model with the one-dimensional transport 
formula^. 

The scattering rate from the initial state (s, fc) is given 
by the T matrix as, 

4 = ^ E Ss',±s, (II) 

sk s'k' 

where S is the system area, and describes the spin- 
conserving scattering time and t~ the spin-flipping scat¬ 
tering time. The momentum relaxation rate and the spin 
relaxation rate are obtained by 
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FIG. 5. Relative spin-relaxation rate Fs/Fp as a function of 
hole density of Uh, calculated for Rs and VmoSs in ih® mono- 
layer M 0 S 2 , VwSa in the monolayer WS 2 . and FmoSs in the 
bilayer M 0 S 2 . Solid lines just represent the linear dependence 
to Uh- 

respectively, where (• • • )gp represents the average on the 
Fermi surface. 

First we present the analyses for Vs and I 410 S 3 in the 
monolayer M 0 S 2 , and also for Vwsa in the monolayer 
WS 2 . Fig.Oplots the relative spin-relaxation rate Fg/Fp 
averaged on the Fermi surface, as a function of the corre¬ 
sponding hole density nh- Here Fp is found to be always 
about ma'^ times a few eV not strongly depending on 
Uh, so it is reasonable to normalize the spin relaxation 
rate by Fp. For HmoSs in the monolayer M 0 S 2 , we see 
that the Fg/Fp is roughly proportional to nh in the low 
density regime, and this is consistent with the qualitative 
estimation Eq. m- The absolute value is of the order of 
10“^ at Uh ^ 10^^ cm“^, which is also in a good agree¬ 
ment with the estimation. This is somewhat surprising 
because Eq. m was derived in the perturbational ap¬ 
proach which is valid in the weak potential limit. 

VwSa in the monolayer WS 2 exhibits a similar behav¬ 
ior, while the absolute value is greater than in M 0 S 2 
about by factor of 10. WS 2 has a larger spin-orbit in¬ 
teraction of AwS 2 ~ 0.211eVi^ In Eq. m, the rela¬ 
tive spin-relaxation rate Fg/Fp is proportional to A^, and 
the order-of-magnitude difference is actually consistent 
with (AwS 2 /^MoS 2 )^ ~ 8.35. For VwSa; we see some non¬ 
monotonic behavior in higher Uh- This may be related 
to the impurity bound states, while the detailed study is 


out of the scope of the present work. 

In a simpler defect Vs, the spin relaxation rate is found 
to be extremely small compared to that of EmoSs j nnd it 
is also deviating from a linear function of nt- In the 
low-energy bands, the wave amplitude of the electronic 
state mainly resides on the transition-metal atoms, and 
thus a single chalcogen vacancy is expected to have rela¬ 
tively small effect compared to EmoSs- seems not 

enough to fully explain the difference of the order of a few 
magnitudes, and we presume some other factor, e.g., can¬ 
cellation of the leading-order matrix elements, may also 
contribute to the suppression of the intervalley process. 

Finally, we study the spin relaxation rate in the M 0 S 2 
bilayer with a EmoSs defect. Unlike monolayer, the en¬ 
ergy bands in the bilayer TMD are all spin degenerate due 
to the coexistence of time reversal symmetry and spatial 
inversion symmetryii However, each eigenstate cannot be 
expressed as a single spin state because the spin-orbit 
interaction strongly hybridizes the different spin compo¬ 
nents in the multiple orbitals. To define the spin relax¬ 
ation rate using Eq. m in this situation, we reconstruct 
the degenerate Bloch states at each Bloch momentum 
into approximate ”spin-up” and ’’spin-down” states so 
as to maximize the expectation value of <Tz . The relative 
amplitude of the spin-down components in the ” spin-up” 
state in bilayer M 0 S 2 is shown in Fig.[4l where we still see 
a strong hybridization about a few of 0.1. In the mono- 
layer, in contrast, the spin hybridization is extremely 
weak as we have already seen, and this is because the 
selection rule imposed by the mirror reflection symmetry 
allows the spin mixing only between the low-energy even 
bands and the high-energy odd bands. In bilayer, the 
mirror reflection symmetry is absent and there the spin- 
orbit interaction directly couples the opposite spin states 
in the same low-energy bands, giving a strong spin hy¬ 
bridization. The spin relaxation rate Fg/Fp calculated for 
bilayer M 0 S 2 is plotted in Fig. [5l and it is much greater 
than that of monolayer by the factor of 10^ reflecting the 
strong spin mixing effect. 

To conclude, we studied the electronic spin relaxation 
in hole-doped TMD monolayer and bilayer in the pres¬ 
ence of realistic atomic defects. By analyzing the band 
structure and spin-orbit interaction, we qualitatively de¬ 
scribe the spin relaxation mechanism and actually esti¬ 
mated the spin-relaxation rate for several specific cases in 
the numerical calculations. In the monolayer TMD, the 
inter-band spin hybridization is suppressed by the parity- 
selection rule under the mirror symmetry, resulting in a 
relatively smaller spin relaxation rate. The bilayer TMD 
has a much greater spin relaxation rate in contrast, be¬ 
cause of the strong spin hybridization in the absence of 
the mirror symmetry. 

This work was supported by Grants-in-Aid for Scien- 
tihc research (Grants No. 25107005). 
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